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Abstract 

We study the Casimir problem for a fermion coupled to a static background 
field in one space dimension. We examine the relationship between interactions 
and boundary conditions for the Dirac field. In the limit that the background 
becomes concentrated at a point (a "Dirac spike") and couples strongly, it 
implements a confining boundary condition. We compute the Casimir energy 
for a masslike background and show that it is finite for a stepwise continuous 
background field. However the total Casimir energy diverges for the Dirac spike. 
The divergence cannot be removed by standard renormalization methods. We 
compute the Casimir energy density of configurations where the background 
field consists of one or two sharp spikes and show that the energy density is 
finite except at the spikes. Finally we define and compute an interaction energy 
density and the force between two Dirac spikes as a function of the strength 
and separation of the spikes. 



1 Introduction 

When a quantum field is coupled to matter, the spectrum of its quantum fluctuations 
changes. The resulting energies, forces, and pressures are known in general as Casimir 
effects and have a variety of experimental consequences, |2] In many applications 
the details of the coupling to materials can be idealized by boundary conditions, and 
the results depend on the geometry and the boundary condition alone. The Casimir 
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force between grounded metal plates has now been measured quite accurately and 
agrees with Casimir's original prediction [2] . 

In Refs. HI new methods were developed to examine Casimir effects without 
resorting to boundary conditions ab initio. Instead the fluctuating fields are coupled 
to a non-dynamical background field, the renormalized effective energy is calculated 
using standard methods of quantum field theory and renormalization theory, and ex- 
amined in the limit that the background becomes sharp and strong enough to enforce 
a boundary condition on the fluctuating fields [5]. In general the renormalized effec- 
tive energy diverges in this limit, indicating that the physical situation depends in 
detail on the coupling between the field and the matter. In the examples studied in 
Refs. 12111] forces between rigid objects and energy densities away from boundaries re- 
main finite and agree with the boundary condition idealization. However the Casimir 
pressure on an isolated surface was found to diverge in the boundary condition limit, 
signalling that it depends in detail on the properties of the material that provide the 
physical ultraviolet cutoff. ^ 

In this paper we apply the methods of Refs. [3 Ej to fluctuating Dirac fermions 
in one dimension. Fermionic Casimir effects are interesting in their own right, for 
example in the bag model of hadrons. [71 |H] Boundary conditions are in general more 
disruptive to Dirac fields than boson fields because the equations of motion are first 
order, so we anticipate that fermion Casimir effects will provide an instructive labo- 
ratory to study the difference between the approach of Refs. [3 Ej and the boundary 
condition method. 

Our object is to treat the boundary condition as the limiting form of the cou- 
pling to a smooth, finite background potential, £int = \il'V{x)xl'. We concentrate on 
potentials, V{x), that couple like the mass, because such backgrounds lead to a con- 
fining (bag) boundary condition when they become strong. [7J As in the scalar case, 
the boundary condition limit is achieved by letting V{x) become concentrated at one 
point (the "sharp" limit) - we refer to this standard configuration as a "Dirac spike" 
- and then letting A — > cxo (the "strong" limit). In the Dirac case the renormalized 
vacuum fluctuation energy diverges already in the sharp limit. In the boson case it 
was finite for a sharp background in one dimension and diverged only in the limit 
X ^ oo. Our results on the energy density and the force between "Dirac spikes" 
are more positive. We find that the renormalized vacuum fluctuation energy density 
away from the spikes is finite even in the strong limit, and that the force between 
Dirac spikes is also finite. 

In Section |21 we discuss the behavior of the Dirac equation in the presence of an 
interaction concentrated at a single point. We show that it is most convenient to 
represent the interaction by a transfer matrix which relates the Dirac wavefunction 
on the immediate left of a singularity to the wavefunction on the immediate right. In 
this way we are able to define the Dirac spike which is the subject of much of the rest 
of our work and which implements a bag boundary condition in the limit X —>■ oo. 

In Section |21 we briefly review the method developed in Ref. jU] to compute the 
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total renormalized Casimir energy in a background field. We emphasize that the di- 
vergences that appear in the case of a smooth background are cancelled by allowed 
counterterms already present in the continuum Lagrangian. We then apply this for- 
malism to compute the energy associated with a single Dirac spike, and find that it 
is ill defined. 

In Section m we consider energy densities and forces between Dirac spikes. First 
we introduce a method based on the Greens function jl], which allows us to compute 
the energy density outside the region of interaction. We then apply this method to 
compute the energy densities of configurations with one Dirac spike at the origin, and 
two spikes separated by a distance d, both of which we find to be finite even as A — > cxd. 
We then introduce an interaction energy density, defined as the energy density of the 
two spike configuration minus twice the energy density of a single spike, placed to 
cancel the local divergences in the energy density, and conclude by integrating the 
interaction energy density to compute the total energy. From that we derive the force 
between the two Dirac spikes in one dimension. 

The problem of vacuum fluctuations of the Dirac field in one dimension has been 
studied previously jTUl E]. In Ref. [121113 boundary conditions were imposed ab ini- 
tio, and zeta function regularization was used to control divergences. In Ref. |14|. 
zeta function regularization was used to derive an expression for the total energy of a 
configuration of two parallel plates, also using boundary conditions. In ^^1; expres- 
sions for zero-point energies and energy densities were developed using a zeta func- 
tion regularization framework, but the physical potential implementing the boundary 
conditions was not taken into account. In this paper we implement the boundary 
conditions as the limit of a finite potential. In this way we can control the subtleties 
of the calculation at every step enabling us to confidently distinguish what is finite 
and cutoff independent and what is divergent in the boundary condition limit. 

This work is related to Ref. , which applied the methods of Ref. jH] to fermion 
fluctuations in a smooth background in one dimension, but did not consider either the 
energy density or the limit of sharp and strong backgrounds necessary to implement 
boundary conditions. In this paper, we will adopt the renormalization conventions of 
Ref. and refer the reader there for further discussion. 

2 The Dirac Equation in a Singular Background 

The equation of motion for the fluctuations of a Dirac field in a non-dynamical back- 
ground field, V , in one dimension reads 



Here, a = 7o7i, = 70 are two-by-two matrices, and ui is the energy eigenvalue of 
the time-independent solution ip{x) given the potential V{x). We work in a chiral 
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basis, a = as, P = (Xi, 75 = 7071 = (13. Our spinor normalization for right- {k > 0) 
and leftgoing {k < 0) waves (when y = 0) is determined by 



iPoik, x) = t;(fc)e'^" = ( ^/^±? ) e"'^ (2.2) 

yv LO — k J 

with k = — m^. Note that V"(x) is matrix- valued. In its most general form, 

V{x) = Vo{x)I + Vi{x)a + V2{x)p + V3{x)ap , (2.3) 

where Vq, Vi, and V2 are real and V3 is imaginary. For simplicity we exclude back- 
grounds proportional to a(3. This eliminates interactions of the form Ci^t = '^'l^^si^)'^ ■ 
In addition, Vi can be gauged away. To show this, assume that ip{x) solves eq. ()2.H) 
when Vi{x) = 0. Then it is easy to show that 

Xix) = A{x)iIj{x) = e-'^^ ^i("')'^">(x) (2.4) 

solves eq. ()2.H) for any V\(x). Since V\ only affects the phase of the solutions (and 
is independent of energy), the density of states as a function of x and lo is indepen- 
dent of Vi, and we conclude that the energy density and hence the total energy are 
independent of V\. We can therefore, without further loss of generality, consider only 
potentials that can be written as a linear combination of the identity and /3. To 
demonstrate the physical significance of the two remaining terms, reorder the terms 
in eq.(EIID, 

d 



-ia— + (3(m + M(x)) 
ax 



ipix) = {io -(j){x))'^{x). (2.5) 



where we have replaced Vo{x) by (f){x) and ^(x) by M{x) 

The component proportional to (3 in effect acts like an extra mass term and the 
component proportional to the identity acts like an electrostatic potential. The mass 
term affects positive and negative energy eigenstates identically, both are attracted 
or both repelled. The electrostatic potential, on the other hand, treats positive and 
negative energy eigenstates in an opposite manner, attracting one and repelling the 
other. A significant portion of the later discussion will be focused on the pure mass- 
like potential. 

As explained in the Introduction, we intend to implement a boundary condition 
as the limit of a finite potential. We must first investigate how this can be done for 
the Dirac equation in a consistent way. For simplicity, whenever the explicit form of 
the finite potential is needed, we will use a rectangular barrier. As the barrier width 
a goes to zero while its area remains constant, the potential can be replaced by a 
transfer matrix relating at the left of the barrier to at the right. We call this 
object a "Dirac spike". If we then let the area under the barrier go to infinity, the 
domains to the left and right of the spike decouple and obeys a flux conserving 
boundary condition at the spike. Depending on the matrix nature of the potential, 
the Dirac spike may have up to two bound states. 
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In the scalar case the representation of a pointhke potential is simple: a delta 
function potential serves the purpose. Since the Dirac equation is first order, however, 
a delta function potential leads to contradictions. For example, consider a masslike 
potential V{x) = (3M5{x). Written out in components the Dirac equation becomes 

- 2^ + (m + M5(x))^2 = i^i^i (2.6) 
dx 

+ (m + M5(x))V'i = (2.7) 

dx 

The terms involving a delta-function are only well defined if ip is continuous at x = 0. 
However the first equation implies a jump in ipi for continuous 1^2, and the second 
requires a jump in ip2 for continuous ipi. Thus the equations are not consistent. 

Anticipating that a sharp potential will generate a discontinuity in ip at the point 
of interaction, we study a sharp potential as the limit of a finite one. Consider a 
potential of the form 

V{x,a) = Tb{x,a), (2.8) 

where b{x, a) is defined to be 1/a for < x < a and otherwise, and F is some 
matrix- valued constant. It is easy to show that the solution to the Dirac equation in 
a potential of the form of eq. ()2.8|1 is a continuous function of x. 

Since the Dirac equation is first order, iIj{x) is completely determined from its value 
at any point. In particular, given ip{0), we can compute ip{a), ip{a) = T(T,a,u)ip{0), 
where the two dimensional transfer matrix, T, depends on uj, a, and F. As long as 
a is finite, ip is continuous everywhere. But as we take the limit a 0, keeping F 
fixed, ip is no longer continuous at x = 0. We therefore define the transfer matrix for 
a Dirac spike by T(F, uj) = lima_^o T{r, a, u) with the jump condition 

^(0+)=T(F,a;)^(0-). (2.9) 

where ip{0'^) = lim^.^o± 'ipix). 

The properties of T depend on both the strength and Dirac matrix character of 
F. As advertised we limit our discussion to potentials of the form, 

T = ei + X(3, (2.10) 
where 6 = (pa and A = Ma are constants. For finite a we obtain 

I cos qa + i^^^ sin qa _j^ rn+M ^-^^ \ 
1 smga cos qa — i—^ sm qa j 

where q = sj {lo — (pY — (m + M)^. In the limit a — for any fixed value of ui, we 
obtain the transfer matrix for the general Dirac spike. 



T{X,9) 



cos V0^ - A2 - i^^^ sin _ ^2 _i_^sm _ p 

i-^=^sin V^2_p cos V0^ - A2 + i^^S^ sin V0^ - 

(2.12)' 
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This is the generahzation of the delta function potential to the Dirac equation (modulo 
our exclusion of pseudoscalar couplings). Note that in the limit a — T becomes 
independent of u and can therefore be regarded as a background "potential" in the 
Dirac equation. 

Two special cases are of interest. The masslike case, where 6 = 0, will play the 
principal part in the sections that follow, 

Tm{X) = ( -.sinh A\ 
\t smh A cosh Ay 

For completeness, we will also compute the electrostatic case, where A = 0, 

Te(^)=('q' j.)=e-^^^^ (2.14) 

These two potentials may or may not have bound states. For the pure mass-like 
case, there are positive and negative energy bound states when A < 0. The energy 
eigenvalue obeys 

^!E<^^ = _tanhA, (2.15) 
m m 

whereas in the electrostatic case, there is one bound state in either the positive or 
negative energy spectrum regardless of the sign of 6, 

^^^^^^ = -i.ne (2.16) 

Equations ()2.15p and ()2.16p are qualitatively different. Whatever the sign of the 
electrostatic interaction, it is attractive either for particles (positive energy states) 
or antiparticles (negative energy states). The electrostatic interaction violates charge 
conjugation invariance so the spectrum need not be symmetric, however it is periodic 
m 9 ^ 9 + 2ti. The masslike spike has either no bound states when it is repulsive, or 
exactly one positive energy and one negative energy bound state when it is attrac- 
tive. The energy spectrum is symmetric since a masslike interaction preserves charge 
conjugation invariance. 

The probability current passing through a masslike spike goes to zero as the 
strength of the spike goes to infinity. Effectively the spike becomes a wall, split- 
ting the line into two independent half lines. This situation is familiar from bag 
models, and not surprisingly the masslike spike reduces to a bag boundary condition 
in this limit. |7. In contrast the electrostatic spike is periodic in the strength 9 and 
therefore does not reach a limit as — oo. While this phenomenon is interesting in 
its own right, we do not pursue it any further here. 

As A — > oo, eq. ()2.13|) becomes 

Tu{\) = e'^(\ 7)+0(e-^) (2.17) 
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For '?/'(0"^) to be finite, it is necessary that 



^i(o-) 



#2(0-) + 0{e-^) 



(2.18) 



from which it follows that 



V'i(o+) 



^V'2(0+) + O(e 



-A 



(2.19) 



Eqs. (I2.18P and (j2.19|) can be written as two independent boundary conditions on the 
half lines separated by a true boundary at x = 0, 



which are standard bag boundary conditions and insure that the probability current, 
j = ipjiip vanishes at a; = 0. 

3 Total Renormalized Vacuum Fluctuation Energy 

In this section we briefly review the method developed in Ref. [HI 0] to compute the 
Casimir energy of a field configuration from scattering data, and adapt it to the case 
of a Dirac fermion in a static background in 1+1 dimensions. Our work is based on the 
prescription of Ref. |16], specialized to a sharply peaked background. We first review 
the fundamentals of the method, including how the phase shifts are related to the 
density of states. We then show how the resulting highly oscillating energy integrals 
can be evaluated using contour integration, significantly simplifying the numerical 
computations. 

We first confront the case of a Dirac spike directly and find that it generates a 
divergent renormalized vacuum fluctuation energy. We then compute the total energy 
of a potential of the form of eq. ()2.8|1 . only afterwards taking the potential to be sharp 
and strong. Using this approach, the total energy is well defined for finite width, a, 
of the potential, but diverges as a — > with the integral of the potential fixed. 

In one dimension the S matrix can be written in terms of the transmission and 
reflection coefficients. 



The transmission coefficients for right and leftgoing waves must be identical by time 
reversal symmetry, but the two reflection coefficients can differ. The S matrix is 
constrained to be unitary, S'^S = I. In addition, as k becomes large, any finite 
potential becomes negligible, and we expect the 5* matrix to become the identity in 
this limit. Of course we are considering singular potentials, so this property must 



(1 + Z7i)^(0+) 



0, 







(2.20) 
(2.21) 




(3.1) 
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be studied carefully. In a basis of eigencliannels, the S matrix is diagonal, and its 
eigenvalues are complex exponentials of the phase shifts, 

/g2i5i(fe) Q \ 

•S" = ( Q g2i<52(fc) ) • (3-2) 

The change in the density of states, Ap{k), due to the interaction is related to the 
trace-log of S, A{k) = ^Trln5(A;) = 6i{k) + 62{k) f^, 

A /,N /,N n/,N I dA(k) 1 rflndetS' ,^ ^, 

Following Refs. [S] and JH] we write the renormalized Casimir energy for any 
non-singular background, V{x), as 

POO 

AEcas = - y"(c^i - m) - / rffc(cu(A;) - m)Ap{k) + TpD, (3.4) 
J Jo 

Note the negative sign and factor of two relative to the Casimir energy of a real scalar 
field. The sum is over possible bound states, and 

ApW = iAs(,)m.(A„)_Aa.„)), (3.5) 

where A^^\k) is the first Born approximation to A{k). The subtraction of the first 
Born approximation is compensated by adding back the corresponding Feynman di- 
agrams, which are then combined with the counterterms to give the renormalized 
Feynman diagram contribution, Ffd- We work in the "no tadpole" renormalization 
scheme where the counterterms exactly cancel the local, divergent Feynman diagrams. 
In this renormalization scheme FpD = 0. For more discussion of the derivation of 
eq. (|3.4p we refer the reader to Ref. [TB] . 

The formalism of Ref. [Sj and eq. ()3.4|) gives a finite renormalized energy for any 
finite potential. However, the total energy of a Dirac spike is infinite. This is easy 
to see by computing the associated shifted density of states and and examining its 
behavior at large k. 

The S matrix can easily be obtained from the transfer matrix, eq. ()2.13j) . The 
result is 

S=(l = V..,...., ( . ,i„ , -'"f^") . (3.6) 



R2 T J k cosh \ + im sinh A \-iuj sinh A k 

It is already apparent that the Dirac spike is a badly behaved. We expect 5* ^ / as 
/c — * 00, but instead 

^ f sechA — ztanhA^ 
\ —i tanh A sechA 
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Furthermore, subtraction of the first Born approximation, ie. of the terms hnear in 
A, does not suffice to render the k integral convergent. Specifically 

. , , , 1 , / , N 1 ^ sinh A cosh A , „, 

Ap(A;) = -— A(A;) = 5 5—. (3.7) 

Ti dk IT k^ cosh A + sinh A 

so the integrand in eq. (j3.4j) . {uj{k) — m)Ap{k), is proportional to ^^^^^-^ at large 
k, and the k integration diverges logarithmically. It would be necessary to subtract 
a term proportional to tanhA in order to render the integral finite, and there is no 
justification in conventional field theory for a subtraction that includes all orders in 
A. We conclude that the Casimir energy diverges for a Dirac spike, even at finite A. 

Having encountered a problem confronting the Dirac spike directly, we back off and 
instead compute the total renormalized Casimir energy of a finite mass barrier, and 
study the difficulties that develop as the potential becomes sharp. We first compute 
the shift in the density of states by calculating A{k) = ^Indet S. The S matrix for 
this configuration is 

^ _ 1 f kq e~'*^° —iMuj sin qa\ ,^ g, 

~ kq cos qa + i{mM - k"^) sin qa \-iMuj sin qa e'^*'"' kq e~''"' ) ^ ' ' 

From this we conclude that 

^ A{k) ^ ^ \ dtS aM {k^m + k'^M cos'^ qa — rn^M sin^ qa) 
dk 2i dk k'^q^ + M^cj^ sin^ qa 

sin qa cos qa{k'^ + niM + 2m?) 



k'^q^ + M'^uj'^qsiY? qa ^ ^ 

Taking M — cx), a — >■ 0, with aM = A fixed, sends qa ^ iX, and eq. ()3.9p becomes 
eq. ()3.7|) . as it must. The first Born approximation is given by^ 

A^'\k) = - — {m+-). (3.10) 
We compute the Casimir energy from eq. ()3.4j) 

EcUm,M,a) = -- {u: - m)—{A{k) - A^^\k)) dk 

J 

/-"^ dk(uj-m) f k^q 

cos2ga 



^ Jo k^Q^ + k'^qi^'^M'^ sin qa \ 2 
+ k^li.^ + 2mM + 2m^) - q sin^ qa{uj'^M{m + ^) + m'^k'^) 

— A;^ sin ga cos ga(/c^ + mM + 2m^)^ (3-11) 

It is easy to see that Ecas{fn, M, a) is finite for finite M, as it must be, but it diverges 
as AM In ^ as a ^ with X = Ma fixed. 

m 

^The first Born approximation in the Dirac theory is equal to the tadpole plus the local part of 
the self-energy Feynman Diagram. The explicit form can be derived by writing the Dirac equation 
as a second order equation and applying the general formalism. |16| 
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4 Greens Functions and the Energy Density 

The phase shift formahsm does not allow us to compute energy densities. In this 
section we extend the formalism of Ref. [4^ based on the Greens function to enable 
us to study energy densities. We relate the trace of the Greens function to the local 
density of states and integrate over k to obtain the local energy density. We then 
construct the Greens function from solutions to the Dirac equation. We compute 
the free Greens function and the interacting Greens function for one and two spike 
potentials, and finally use the result to compute the force between two Dirac spikes. 
In one dimension, the Greens function S{x, x',u) is a 2 x 2 matrix, and obeys 



-00 — ia- — h V(x))S(x, x , uo) = 6(x — x) (4.1) 
ax 



We can write 5* as 



where 



^(x,x-,c.) = V ^"("^^"^"^^ (4.2) 
^-^ ujr, — ijJ — le 



{-icy-y- + V{x))'ljjn = tOnl/jn (4.3) 

dx 



and 



J2Mx)^Ux') = (I ^^S{x-x'). (4.4) 
From eq. (|4.2|) we define a local density of states, 

p{x,uj) = - Im Tr S{x,x,u) = S{un - u;)V^^(x)?/;„(x) (4.5) 

TT ^ — ' 

n 

We can use the local density of states to compute the local energy density (or indeed 
any other local density functional of the Dirac field). To compute the energy density 
we integrate the density of states weighted by the energy. However, we must remember 
that the ground state energy for a fermion is —u, so 

e(x) = / da; a; Im Tr S{x,x,uj) = / dk k Im Tr S{x, x,uj{k)) (4.6) 

Jo ^ io 

where S{x,x',uj) = S{x,x',uj) — So{x,x' ,uj), where 5*0 is the free Greens function. 
The free Greens function 5*0 is given by 

So{x, x', t^) = y - x)v{-k)v\k)(3 e''=("'-") + 

e{x - x')v{k)v\-k)l3 e*fc("-"') } . (4.7) 
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where v{k) is the free Dirac spinor defined by eq. ()2.2|) . For a configuration with one 
spike, the Greens function is 

= S„(..x'..) + ^^-^^^^^ .(k)Ak)0 e"<-') (4.8) 

and for a configuration witli two spikes located at a: = ±(i/2 it becomes 

S2ix,x',uj) = So(x,x',uj)-—T2(uj)smh^\e^''"^x 



2 

— T2(cj)(cosh A + — sinh A) sinh A e''"^ x 
k^ k 



v{k)v\k)f3 e^'^(^+^') + v{-k)v\-k)f3 g-^'^^^+^'M , (4.9) 



between the spikes and 

S2ix, x', u) = Soix, x', u) + ^^^^^vik)v\k)f3 e*'=("+"') x 

rh 

{k sinh 2 A cos kd + 2m sinh^ A sin kd) (4-10) 

outside the spikes. Here T2{uj) is the transmission coefficient for the entire two-spike 
potential, given by 

k'^ 

T2{uj) = 2 — (4-11) 

(k cosh A + im sinh A)^ + sinh A e'^^'"^ 

We combine these expressions for the Greens functions with the definition of the 
energy density, eq. ()4.6|) to obtain the energy density for the one and two Dirac spike 
configurations. Note that no subtractions are needed away from the spikes. The 
energy densities are finite, and anyway, the subtractions are local in the background 
field, and therefore vanish except at the spikes. 

For the single spike configuration, it follows from eq. ()4.8|) and the symmetry 
under x — >• —x that 

Tr{S,{x,x,u;)-So{x,x,u;)} = ^ f^^^^T+L^^nh X ^^''^^ 

for all X ^ 0. Using eq. ()4.6|) . we can now write down an expression for the energy 
density at x 7^ 0, 

1 ^ mu; sinh A e^*''!''! 



The Casimir Effect for Fermions in One Dinaension 



12 




Figure 1: Energy density for configurations with a Dirac spike of strength A = 0.1 
(sohd), A = 0.5 (dashed), and as A — >• oo (dotted). The energy density converges 
uniformly as A ^ cxd. 

For X 0, the integrand of eq. ()4.13j) falls off exponentially in the upper half /c-plane, 
but oscillates rapidly on the real axis. It is therefore more readily evaluated using 
the method of contour integration. The integrand has a branch cut on the positive 
imaginary axis, but no poles in the upper half plane because there are no bound 
states, thus 

m2 V^^^ sinh A e-^-l^l 



^^^^ TT J I (r cosh A + sinh A) ^ 

where ^ = mx. Figure ^ shows a plot of the energy density for several choices of A. 
The divergence at x = is clearly visible for all values of A. 

As a check on this method we compute the integrated change in the density of 
states Ap(/c) by integrating p{x,u) over x. 

/•oo ^ r-oc, m sinh A e^**^!^' 

Ap(fc) = / dx p(x,k) = — Im- — ——dx. (4-15) 

vr J__^ k cosh X + im smh A 

The integrand is even in x, so we evaluate it only for positive x, and then multiply the 
result by two. At the upper limit, we take x — > oo(l + ze), and get no contribution, 
yielding the result 

1 muj sinh A cosh A 

vr A;2 cosh^ X + iv? sinh'' A 

identical to eq. (j3.7|) as expected. 

When we turn to two spikes, we maintain parity invariance by giving the two 
spikes the same strength parameter A, and placing them at x = ±(i/2. The energy 
density is even under x —x. Hence, there are only two distinct regions: the region 



= , . ■ .2. - (4-16) 
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Figure 2: Energy density in units of for a configuration of Dirac spikes of strength 
A = 1 of separation d, with m = 0.1/ d (dotted), m = 1/d (dashed) and m = 10 /d 
(sohd). 



between the spikes, and the region outside them. Written as integrals, the energy 
densities are 



For \x\ < d/2, e2{x,d,X,m) 



For \x\>d/2, e2{x,d,\,m) 



■- I 



mVt'^-m'^{t sinh 2X+2m sinh^ A)e-*'^ cosh 2t\x\ . 
m {t cosh A+msinh \)'^+(t'^—m?) sinh"^ A e-^td 



1 

TT 



Jrr 



(t cosh A+m sinh X)'^ +(t'^ —m?) sint? A 



mf-y/t-^—m'^ sinh 2A cosh fd e 



(t cosh A+m sinh A)2+(t2_m2) sinh^ A e-2*d 

2m-^\/t-'^— m-^ sinh^ Asinhf:d e"-^*!"^! 
(t cosh A+m sinh \)'^+(t'^—m?) sinh"^ A e"^*'' 



+ 

(4.17) 



The integrals converge quickly except at x = ±(i/2, where they diverge. A plot of 
the energy density for several choices of m is shown in Figure El 

To eliminate local divergences which do not affect the force, we consider the pure 
interaction energy density e(x, (i), defined as 



e(x, d) = e2{x, d) — ei(x — d/2) — ei(x + d/2), 



(4.18) 



where ei{x — d/2) is the energy density of a single spike placed at x = d/2 defined 
in eq. ()4.14|) . Because the divergences are local to the points of interaction, this 
subtraction will render the resulting expression finite everywhere. However, there is 
no reason to expect that the energy density should be continuous across the point 
of interaction, and this is indeed not the case. In particular, for m = 0, the energy 
density vanishes in the region outside the spikes, but takes a constant nonzero value 
in the region between them. Figures |3] and El show the interaction energy density in 
various configurations. 
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e(x,d,A,,m)d' 




x/d 



-0.1 [ 

Figure 3: For small m, the interaction energy density is significant, and essentially 
localized between the two spikes. As m grows larger, the overall interaction term gets 
smaller, and energy density "leaks out" to the region outside the spikes. The plot 
shows m = (solid), m = 0.1/d (dashed), m = 0.3/d (finely dashed), and m = 1.0/d 
(dotted), with A = 1 fixed. 

e(x,d,X,m)d^ 



0.005 



-0 



-0.005 



-0.01 



x/d 



Figure 4: Interaction energy density as a function of A. The plot shows A = 0.4 (solid), 
A = 1.0 (dashed) and A ^ oo (dotted), with m = 1 fixed. Note the discontinuities at 
X = ±d/2. 
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Finally we compute the force between two spikes. Although the total Casimir 
energy diverges, the divergences are localized to the points of interaction and are 
identical to the divergences associated with a single spike. Therefore we define an 
"interaction Casimir energy" by subtracting away twice the contribution from a single 
spike from the two spike Casimir energy. The subtracted quantity is independent of d 
and therefore does not contribute to the force. To simplify the computation we place 
the centers of the two subtracted spikes at x = ±(i/2 respectively, and perform the 
subtraction within the integrand of the expression for the energy. As demonstrated 
in Sectional this renders the energy density finite everywhere. 

We define the total interaction energy E as the integral over x of e(a;, d), defined 
by eq. 



-2td ^ 



/oo 1 /"^ 

dx e{x,d) = — dt 2^/t'^ — w? sinh^ A e~ 
-oo ^ Jm 

{w?{td + 1) — t^d) cosh A — [mdif' — w?) — mi) sinh A 



(t cosh A + m sinh A)^ + (t^ — m?) (t cosh A + m sinh A)e 



2td 



(4.19) 



Since this integral is absolutely convergent, we can compute the force F{X,m,d) = 
—dEiat/dd by moving the derivative operator under the integral sign. The resulting 
expression in not very enlightening, and the t-integration must be performed numer- 
ically, however in certain interesting limits, the formula simplifies greatly. 

In the limit A — oo, corresponding to impenetrable spikes, or equivalently, bag 
boundary conditions, the total interaction energy, eq. ()4.19|) reduces to reduces to 

hm Eint(A, m, d) = E^^^{ni, d) (4.20) 

A-^oo 



2 ^- rn-dit^-m^) 

- / - m2- ^r^dt 4.21 



and the force simplifies as well, 

lim F(A, m, d) = F^^g{m, d) (4.22) 

A^oo 

2 -e'^^'^Um'^ - t^) + 2t3rf - 2mt{md + 1)1 -it- mf , 

= / Vt^—m^ r ; , dt 

71 Jm [{t -m) + {t + m)e2*^]2 

(4.23) 

These integrals can be done analytically for m = 0, 

^bag(0,rf) = (4.24) 
^bag(0,ci) = (4.25) 

an attractive inverse-square dependent force law. 



The Casimir Effect for Fermions in One Dinaension 



16 
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Figure 5: Force as a function of separation between the spikes. We write F 
— ^^f{X,md), and the plot shows / for A = 0.5 (sohd), A = 1 (dashed), and A 
(dotted). 



oo 



In general, the force can be parameterized by 

F{X,m,d) = -^f{X,md), (4.26) 

where / is a dimensionless function of two dimensionless parameters and describes 
the departure from the force in the massless case. Figure El shows a plot of f{X,md), 
and Figure El shows /(A,0). f{X,md) approaches a limit exponentially as A ^ oo, 
/(A, 0) ~ f{oo, 0) + C e~^. As a consequence, all derivatives of / go to zero as A ^ oo. 
This is visible in Figure IHl 



5 Discussion 

In this paper we have applied the methods of Ref. |3] to fermion fields in one dimen- 
sion. The background fields were taken to be sharp spikes, imposing jump conditions 
on the dynamical fermion field. Using the method of computing the density of states 
directly from the S matrix, we showed that the total energy of such a configuration 
diverges and thus cannot be computed directly using the jump condition approach. 
The details of the interaction cannot be ignored, and the divergences introduced by 
a sharp potential in a fermion field are even worse than the ones that appear when 
the dynamical field is a scalar. 

Nonetheless, energy densities and forces are still well defined. Using a formalism 
based on the Greens function, we showed how to compute these quantities even when 
the background field is highly singular, and derived expressions for the energy density 
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f(1/^,0) 




2 4 6 8 10 



Figure 6: For m = 0, the force takes the form F = ~2iipf(^^ dimensional reasons. 
This plot shows f{l/^) as a function of ^ = 1/A. 



in backgrounds of one and two spikes. The work culminated in a computation of the 
force between two spikes in the presence of a fermion field. For a massless field and 
spikes implementing bag boundary conditions, dimensional considerations require the 
force to be an inverse square force as a function of the distance between the spikes. 
We explicitly verified this, and computed the constant of proportionality. 
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